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^ ' Abstract 

, We review recent progress in the study of varying constants and 

attempts to explain the observed values of the fundamental physical 
■ constants. We describe the variation of G in Newtonian and relativis- 

. tic scalar-tensor gravity theories. We highlight the behaviour of the 

I isotropic Friedmann solutions and consider some striking features of 

O^' primordial black hole formation and evaporation if G varies. We dis- 

5-H . cuss attempts to explain the values of the constants and show how 

we can incorporate the simultaneous variations of several 'constants' 
exactly by using higher-dimensional unified theories. Finally, we de- 
^ ■ scribe some new observational limits on possible space or time varia- 

H ! tions of the fine structure constant. 



1 Introduction 

In this overview of some aspects of varying constants we will begin by con- 
sidering the time variation of the gravitational 'constant' G in Newtonian 
and relativistic theories of gravity. Although the Newtonian situation is usu- 
ally ignored, it provides a number of instructive parallels and contrasts with 
the more complex situation that prevails in scalar-tensor generalisations of 
general relativity. We will focus upon the behaviour of the cosmological solu- 
tions in these theories and provide prescriptions for generating the isotropic 
Friedmann solutions to any version of a scalar-tensor gravity theory. 
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Next, we shall highlight the unusual situation that seems to be created if a 
primordial black hole forms in the very early stages of a universe in which G is 
changing with time. Then we shall go on to consider some of the speculative 
ideas that have been put forward to explain the values of the constants 
of Nature. We shall discuss the problem of the simultaneous variation of 
several 'constants' and describe how this situation can be modelled using 
simple scaling invariances of physics or by exploiting the structure of unified 
higher-dimensional theories of the fundamental interactions. One of the most 
interesting quantities to appear in these discussions is a, the fine structure 
constant. In the final section, we describe some new observational limits on 
any possible space or time variations in the fine structure constant that can 
be deduced from spectroscopic observation of molecular and atomic hydrogen 
absorption lines from the gas around radio-loud quasars. 

2 Some Background to Varying G 

The study of gravitation theories in which Newton's gravitational constant 
varies in space and time has many motivations. It began in 1935 with the 
proposal by Milne of a theory of gravitation with two time standards (one 
for gravitational processes, the other for atomic processes) in which the mass 
within the particle horizon, M/j oc c'C'^t, remains constant with respect to t- 
time, led to the prediction that G oc t in this time. The idea became of wider 
interest in 1937 with the 'Large Numbers Hypothesis' of Dirac (1937a,b 1938), 
that the ubiquity of certain large dimensionless numbers, with magnitudes 
0(10^^), which were known to arise in combinations of physical constants and 
cosmological quantities (Weyl 1919, Zwicky 1939, Eddington 1923) was not 
a coincidence but a consequence of an underlying relationship between them 
(Barrow and Tipler 1986, Barrow 1990a). This relationship required a linear 
time variation to occur in the combination e^G'^uiM (where e is the electron 
charge, ttin the proton mass, and G the Newtonian gravitation constant) and 
Dirac proposed that it was carried by G oc (Chandrasekhar 1937, Kothari 
1938) This led to a range of new geological and palaeontological arguments 
being brought to bear on gravitation theories and cosmological models (Jor- 
dan 1938, 1952, Teller 1948, Dicke 1957, 1964, Gamow 1967a,b). Brans and 
Dicke (1961) refined the scalar-tensor theories of gravity first formulated by 
Jordan and, motivated by apparent discrepancies between observations and 
the weak-field predictions of general relativity in the solar system, proposed 
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a generalization of general relativity that became known as Brans-Dicke the- 
ory. As the solar system and binary pulsar observations have come into 
close accord with the predictions of general relativity so the scope for a the- 
ory of Brans-Dicke type to make a significant difference to general relativity 
in other contexts, notably the cosmological, has been squeezed into the very 
early universe. However, more general theories with varying G exist, in which 
the Brans-Dicke parameter is no longer constant (Barrow 1993a). These the- 
ories possess cosmological solutions which arc compatible with solar-system 
gravitation tests (Hellings 1984, Reasenberg 1983, Shapiro 1990, WiU 1993), 
gravitational lensing (Krauss and White 1992), and the constraints from 
white-dwarf cooling (Vila 1976, Garcia-Berro et al 1995). The crucial role 
that scalar fields may have played in the very early universe has been high- 
lighted by the inflationary universe picture of its evolution. A scalar fleld, 0, 
which acts as the source of the gravitational coupling, G ~ is a possible 
source for inflation and would modify the form of any inflation that occurs 
as a result of the universe containing weakly coupled self-interacting scalar 
fields of particle physics origin. There have been brief periods when exper- 
imental evidence was claimed to exist for a non-Newtonian variation in the 
Newtonian inverse-square law of gravitation at low energies over laboratory 
dimensions (Fischbach et al 1986) and speculations that non-Newtonian grav- 
itational behaviour in the weak- field limit might explain the fiatness of galaxy 
rotation curves (Bekenstein and Meiscls 1980, Milgrom 1983, Bekenstein and 
Milgrom 1984, Bekenstein and Sanders 1994) usually cited as evidence for the 
existence of non-luminous gravitating matter in the Universe. Most recently, 
particle physicists have discovered that space-times with more than four di- 
mensions have special mathematical properties which make them compelling 
arenas for self-consistent, finite, anomaly-free, fully-unified theories of four 
fundamental forces of Nature (Green and Schwarz 1984). Our observation 
of only three large dimensions of space means that some dimensional segre- 
gation must have occurred in the early moments of the universal expansion 
with the result that all but three dimensions of space became static and con- 
fined to very small dimensions ~ 10~^^cm. Any time evolution in the mean 
size of any extra (> 3) space dimensions will be manifested as a time evo- 
lution in the observed three-dimensional coupling constants (Freund 1982, 
Marciano 1984, Kolb, Perry and Walker 1986, Barrow 1987). The effect of 
this dimensional reduction process is to create a scalar-tensor gravity theory 
in which the mean size of the extra dimensional behaves like a scalar field. 
In particular, low-energy bosonic superstring theory bears a close relation to 
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a particular limit of Brans-Dicke theory (see section 4). 

However, despite these interconnections with modern ideas in the cosmol- 
ogy of the early universe, the theoretical investigation of gravity theories with 
time- varying G is still far from complete and, aside from the solar system and 
binary pulsar observations (Will 1993), there are few general observational 
restrictions on scalar-tensor theories which are clear-cut. 

3 Newtonian Varying G 

We shall begin by investigate Newtonian gravity theories with varying G, 
pointing out the relationships that these simple solutions have to the more 
complicated solutions of scalar-tensor gravity theories. In the past there has 
been very httle discussion of the Newtonian case (see Barrow 1996). The 
exceptions are the rediscoveries of Meshcherskii's theorem (1893, 1949): for 
example, by Batyrev (1941,1949), Vinti (1974), Savedoff and Vila (1964), Du- 
val, Gibbons and Horvathy (1991), McVittie (1978) and Lynden BeU (1982). 
These authors all recognised the equivalence of the Newtonian gravitational 
problem with time-varying G to the problem with constant G and varying 
masses. 

Newtonian gravitation is a potential theory that is derived from the axiom 
that the external gravitational potential due to a sphere of mass M be equal 
to that of a point of mass M. This fixes the potential to be equal to 

$(r) = - + Sr' (1) 

r 

where A and B are constants; A = —GM and B = |A, where A is the cos- 
mological constant of Einstein. This argument shows how the cosmological 
constant arise naturally in Newtonian theory, as it does in general relativity. 

Unless otherwise stated, we shall set the cosmological constant term zero 
[B = = A). In section 3.2 we shall discuss how its interpretation differs 
from a p = —p fluid when G is not constant and prove some restricted cosmic 
no hair theorems. 

Consider the Newtonian N-body problem with a time-varying gravita- 
tional 'constant' G{t). If the N bodies have masses m^ and position vectors 
Tj then 
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Now, if we have a solution, rj{t), of these equations with G — Gq independent 
of time, then, 

r.M = (^fe(-^ + c) (3) 
is an exact solution of the equations (1) with 

GW^G„x(J^) (4) 

where b, c and to are constants with to 7^ 0. Thus, given any solution of a grav- 
itational problem (for example the output from a cosmological N-body grav- 
itational clustering simulation) with constant G we can immediately write 
down an exact solution in which G varies inversely with time. For example, 
suppose we take the simplest Newtonian cosmological model with zero total 
energy, when G = Go is constant. Then, the expansion scale factor of the 
universe is 



f oc (5) 

By the theorem we have that 

r(t) = (^)(-^ + cr/^ (6) 

when G{t) varies as 

G{t) = Go X ; to 7^ (7) 

and so r(t) oc t^/^ as t — > oo. 

The result (4) is also useful for modelling small variations in G over short 
timescales. If we expand an arbitrary analytic form for G(t) to first order in 
t then 

G(t) = Go + Got + ....O(t') ^ Go(l - tGo/Go)-^ (8) 
and this has the form (4) . 
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This result, a consequence of the scale invariance of the inverse- square 
law of force, was first found by Meshcherskii (1893). It has often been 
rediscovered and elaborated. Duval, Gibbons and Horvathy (1991) have 
explored its existence in a wider context and displayed similar invariances 
of the non-relativistic time-dependent Schrodinger equation with Coulomb 
potential (see also Barrow and Tipler 1986) which enables solutions with 
time-varying electron charge (e^ oc t) to be generated by transformation of 
known exact solutions with constant values of e. In the next section we shall 
prove a generalization of Meshcherskii 's theorem for cases where the pressure 
is non-zero and the equation of state has perfect fluid form. 

3.1 Newtonian Cosmologies with G(t) oc 

We adopt the standard generalization of Newtonian cosmology (Milne and 
McCrea 1934, Heckmann and Schiicking 1955, 1959) to include matter with 
non-zero pressure and a perfect fluid equation of state. We shall conflne our 
attention to isotropic Newtonian solutions. This is of particular interest for 
the real universe in the recent past but we also know that anisotropic New- 
tonian cosmological models are not well posed in the sense that there are 
insufficient Newtonian field equations to fix the evolution of all the degrees 
of freedom (there are no propagation equations for the shear anisotropies 
(Barrow and Gotz 1989a)) and this incompleteness must be repaired by sup- 
plementing the theory with extra boundary conditions or by importing shear 
propagation equations from a complete relativistic theory, like general rela- 
tivity (Evans 1974, 1978, Shikin 1971, 1972), or by ignoring the evolution of 
the shear anisotropy (Narlikar 1963, Narlikar and Kembhavi 1980, Davidson 
and Evans 1973, 1977). 

Consider a homogeneous and isotropic universe with expansion scale fac- 
tor r{t). The material content of the universe is a perfect fluid with pressure, 
p, and density p, obeying an equation of state (where the velocity of light 
has been set equal to unity) 



with 7 constant, li G — G{t) then the equation of motion for r{t) is 



p = (7 - l)p; < 7 < 2, 



(9) 




4:TTG{t){p + 3p)r 
3 



(10) 



The mass conservation equation is 
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Hence, we have 



p + 3-(p + p) = 0. (11) 
r 



r 

p = r > 0, constant. (12) 



We shall initially be interested in power-law variations of G{t) of the form 



G{t) = Go Cj) (13) 



so we have 



r = -Ar'^r'-^^^ (14) 
where A is a constant defined by 

_ 47rGotg(37 - 2)r 

3 ^ ^ 

so the sign of A is determined by the sign of 37 — 2, as in isotropic general 
relativistic cosmologies. Hence, accelerating universes (r > 0) arise when 
37 > 2 regardless of whether G varies or not. However, these accelerating 
universes need not solve the horizon and flatness problems in the way that 
conventional inflationary universes do; that depends upon the value of n. 

A generalization of Meshcherskii's theorem can be proved for the case 
with p — {'J — l)p. If r{t) is a solution with G — Gq constant, then (Barrow 
1996), 

r(t) = (^)r(-^ + c) (16) 
with 6, c, and to 7^ 0, constants, is an exact solution of (14) with 

G{t)^GoX f-^V '^^0 7^0. (17) 



These results provide a Newtonian analogue to the conformal properties of 
relativistic scalar-tensor theories. We can draw a number of general conclu- 
sions from them. As i ^ 00 we have 

r{t) ^t, iic^Q, V7 (18) 
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r{t) ^-r{ — ] , if c = and 7 0. (19) 

In particular, if we take the solutions with constant G = Gq to be the zero- 
curvature Friedmann solutions then, when c 7^ 0, we have 

r{i) oc f"/^^, if 7 ^ (20) 

r{i) oc exp[HQi\, Hq constant, if 7 = 0, (21) 
and the solutions with G{t) oc t^'^~^ at large time have the form 



r{t) oc t(37-2)/37 ^ if ^ _^ Q _^ ^ 

r{t) oc t, if 7 7^ 0, c = (22) 



r{t) oc texp 



t, if 7 = 0,Vc. (23) 



These are particular solutions only, of course, and their properties need not 
be shared by the general solutions for a given value of n or 7. The 7 = 
solution, (23), does not exhibit inflation and is asymptotic to the solution of 
the equation f = 0. This is a result of the very rapid decay of G{t) oc t~^. 

When 7 = 4/3 there is no possible time- variation of G which preserves the 
scaling invariance and for other positive values of 7 the expansion is slower 
than in universes with constant G; power-law inflation does not occur in the 
varying-G solutions when < 7 < 2/3. 

Equation (14) describes motion under a time-dependent force for which 
there need exist no time-independent energy integral. Therefore we cannot 
write down a Friedmann equation for f in the usual way. However, there 
exists a class of particular exact solution with simple power-law form: 

r{t) oc i(2-n)/37. ^ _^ (24) 

^ (2-n)(37-2 + n) ^ (2 - n)(37 - 2 + n) 
^ 12nGo^X{3^ - 2)t2-" 127r72(37 - 2)G{t)t^ ^ ' 

so p > requires that 
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(2-n)(37-2 + n) 
(37 - 2) 



> 0. 



(26) 



Clearly, when n = 0, these solutions reduce to the familiar zero-curvature 
(zero energy) solutions of general relativistic (Newtonian) cosmology with 
constant G. They describe expanding universes so long as n < 2. They 
are particular sohitions because they do not possess the full complement of 
arbitrary constants of integration that specify the general solution. Solutions 
of this sort suggest that there may exist more general solutions that behave at 
early times hke a solution of the form (24) with one value oin — rii for t <ti 
and with another value n — n2 ior t > t\. K 'bouncing' solution would have 
77-1 > 2 and 712 < 2, for example, with ti « to- There are many examples 
of scalar-tensor gravity theories with cosmological models that display this 
early-time behaviour (Barrow 1993b, Barrow and Parsons 1996). 

We shall not go on to discuss the general solutions of the Newtonian 
evolution equation and the circumstances under which they approach these 
particular solutions. Such a analysis can be found in Barrow (1996). 

3.2 Inflationary universe models with p = —p 

Scalar-tensor gravity theories have provided an arena in which to explore 
variants of the inflationary universe theory first proposed by Guth (1981) in 
which inflation is driven by the slow evolution of some weakly-coupled scalar 
field. The scalar field from which the gravitational coupling is derived can in 
principle be the scalar field from which the gravitational coupling is derived 
or it can influence the form of inflation produced by some other explicit 
scalar matter field. A number of studies have been made of the behaviour of 
inflation in scalar-tensor gravity theories (Mathiazhagen and Johri 1984, La 
and Steinhardt 1989, Barrow and Maeda 1990, Steinhardt and Accetta 1990, 
Garcia- BeUido, Linde and Linde 1994, Barrow and Mimoso 1994, Barrow 
1995). The non-linear master equation governing the evolution for r{t) has 
interesting behaviour in the inflationary cases where f > 0. 

The particular power-law solutions (24)- (26) with 2 > 7 > expand when 
n < 2. Although they accelerate with time (r > 0) whenever 87 — 2 < 0, the 
expansion only provides a possible solution of the horizon problem when 



2 - n > 37 > 0. 



(27) 
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So, in the case of radiation (7 = 4/3) the horizon problem can be solved if 
n < —2 in the early stages of the expansion. 

In the most interesting case, when 7 = 0, and p is constant, the perfect- 
fluid matter source mimics the behaviour of a slowly rolling scalar field whose 
evolution is dominated by its self-interaction potential. When 7 = the 
master evolution equation (14) is linear in r 

r = -Xr^'r (28) 

with A < 0. We are interested in determining the asymptotic behaviour of 
this equation as t — > 00 for all values of n in order to determine when there is 
asymptotic approach to the usual de Sitter solution that obtains when n — 0. 
The solutions fall into three classes according to the value of n. For n < 2, 
the solutions asymptote towards the WKB approximation as i ^ 00 

r{t) ~ exp I I ; n< 2 (29) 

where the constant Hubble parameter, Hq, is given by 

g;,_,,_4xGot5(3T-2)r ^3^^ 

o 

which is positive for 87 — 2 < 0. 
For n > 2 the asymptote is 

r{t) ^t;n>2. (31) 

In fact, this is a particular case of a stronger result that does not assume 
that G(t) is a power-law. The asymptote r ~ i results whenever G{t) falls 
fast enough to satisfy (Cesari 1963) 



t I G{t) \dt< 00. (32) 
When n — 2 the asymptote is 



r{t) 



(33) 



a 



A = 



-{l + VT+AA) 
SttGoP 



3t^ 



> 0. 
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Thus wc see that if G{t) falls off faster than t'"^ the solutions of eqn. (34) 
approach those of f = and no inflation occurs. By contrast, if G{t) falls 
off more slowly that t^^, grows (n < 0), or remains unchanged (n = 0), then 
inflationary solutions of the form (29) arise. We notice that in the absence 
of G- variation {n — 0) this solution reduces to the well known de Sitter ex- 
pansion {r{t) oc exp{ifoO) familiar in general relativistic models of inflation 
with a constant vacuum energy density or positive cosmological constant. 
When < n < 2 it produces a form of sub-exponential inflation which is 
familiar from studies of scalar-tensor gravity theories with varying G{t) and 
models of intermediate inflation studied in general relativistic cosmologies 
containing a wide range of scalar fields, (Barrow 1990b, Barrow and Saich 
1990, Barrow and Liddle 1993). If n < there is super-exponential inflation. 

There is a Newtonian 'no hair' theorem in the general case where no 
particular form is assumed for the time-variation of G{t) because we can 
make use of the general asymptotic properties of the evolution equation. We 
have already given this result for cases where G{t) falls off faster than t"'^ 
as t ^ oo in eqns. (31)- (32), and as in eqn. (33). In the case where the 
fall-off is slower than we have a WKB approximation 

r(t) ~ c [G{t)]-'^ exp jo; J [G{t)]'^ dt^ (34) 
as i — > oo, where cu^ — 1 and c is a constant, so long as 

t^G{t) — 5> oo as t — oo. (35) 

Clearly, eqn. (29) gives this asymptote in the special case that G{t) oc t~" 
with n < 2 when we choose the a; = 1 part of the linear combination of 
solutions. 

When G is constant in general relativity and in Newtonian gravitation 
the presence of a perfect fluid with an equation of state p = —p is equivalent 
to the addition of a cosmological constant to the field equations. However, 
in Newtonian gravitation with varying G and in scalar-tensor generalizations 
of general relativity, these two are no longer equivalent. If we had included 
a cosmological constant term in the Newtonian gravitational equation (10), 
it would become 

m = .t^mki^ + ^ (36) 

3 3 
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and we can see that the choice p = —p = constant, using (11), only makes the 
first term on the right-hand side of (36) equivalent to the term oc Ar if G(t) is 
constant. Clearly, the late-time may not be dominated by the cosmo logical 
constant term when G{t) oc and n < 0. We will not investigate the 
behaviour for general 7 here, see Barrow (1996). 

4 Relativistic scalar-tensor theories 

Newtonian gravitation permits us to 'write in' an explicit time variation of 
G without the need to satisfy any further constraint. However, in general 
relativity the geometrical structure of space-time is determined by the sources 
of mass-energy it contains and so there are further constraints to be satisfied. 
Suppose that we take Einstein's equation in the form (c = 1) 

Gl ^8nGT^ (37) 

where G^ and are the Einstein and energy-momentum tensors, as usual, 
but imagine that G = G{t). If we take a covariant divergence -a of this 
equation the left-hand side vanishes because of the Bianchi identities, T^"„ = 
if energy-momentum conservation is assumed to hold, hence dG/dx^ — 
always. In order to introduce a time-variation of G we need to derive the 
space and time variations from some scalar field, i/j which then contributes a 
stress tensor Tl^{ip) to the right-hand side of the gravitational field equations. 

We can express this structure by a choice of lagrangian, hnear in the 
curvature scalar R, that generalizes the Einstein- Hilbert lagrangian of general 
relativity with, 

L = -f{iP)R + Uai^df^i^ + 167rL^ (38) 

where is the lagrangian of the matter fields. The choice of the function 

f{ip) defines the theory. When ■?/' is constant this reduces, after rescaling 
of coordinates, to the Einstein-Hilbert lagrangian of general relativity. (We 
ignore, for simplicity the possibility of including a cosmological 'constant' 
term which can now be a function of the scalar field ip). For historical reasons 
scalar-tensor theories have not been written with a gravitational lagrangian 
of this simple form (Bergmann 1968, Steinhardt and Accetta 1990, Holman 
et al 1991) but have followed the formulation introduced by Brans and Dicke 
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(1961). This can be obtained from (38) by a non-linear transformation of ip 
and /. Define a new scalar field 0, and a new coupling function a;(0) by 



= f{^) (39) 



then (75) becomes 



L = -(t)R + ^^^dac^d^ct) + 167rL^. (40) 



This has the Brans-Dicke form. The Brans-Dicke theory arises as the special 
case 

Brans — Dicke : uj{(f)) — constant; f{i/j) oc -0^. (41) 

However, in general, there exists an infinite number of these theories defined 
by the choice of cj(0). One of the reasons for renewed interest in scalar- 
tensor gravity theories of this form is the relationship that exists between the 
gravitational part of lagrangian (i.e. excluding L^) for Brans-Dicke theory 
and the low-energy effective action for bosonic string theory, which can be 
written as (Callan et al 1985) 

Lsst = exp(-2x)(i? + 4x'"x,a - ^^') (42) 

where x is the dilaton field and H'^ = HabcH"'^'^, where Habc is the totally 
antisymmetric 3-form field. If we identify = exp(— 2%) then Lsst is identical 
to the Brans-Dicke lagrangian, (40), when uo = —1. However, differences 
arise in the couplings of the scalar fields to other forms of matter in the two 
theories. 

The field equations that arise by varying the action associated with 
L in (40) with respect to the metric, Qahi and separately, gives the field 
equations. 



G 



ab 



-Stt o;(0) 

, -'■ab ,0 



\a:b - 9ab^^ 



- -^9ab<i 



(43) 
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[3 + 2u;{<l>)] □(/) = SttT: - J{<t,)<t>^^^ (44) 

where the energy-momentum tensor of the matter sources obeys the conser- 
vation equation 

T.f = 0. (45) 

These field equations reduce to those of general relativity when (and hence 
a;(0)) is constant, in which case the Newtonian gravitational constant is 
defined by G = . An interesting feature of these equations is clear by 
inspection: when the trace of the energy-momentum tensor vanishes (this 
includes vacuum and radiation-dominated solutions as important particular 
cases) any solution of general relativity is a particular solution of the scalar- 
tensor theory with (and hence G) constant. 

These equations are also conformally related to general relativity when 
the trace — Q. If we conformally transform the metric 

Qah ^'"^Qab (46) 

with Q, = 0^/^ and define ^ by 

then the conformally transformed theory is general relativity with a matter 
source consisting of a scalar field ^ with a potential V{^). This conformal 
invariance can be exploited to produce powerful solution-generating proce- 
dures for scalar-tensor theories (Barrow 1993a, Barrow and Mimoso 1994, 
Barrow and Parsons 1996). One can regard Meshcherskii's theorem and its 
generalization proved above in (16)- (17) as Newtonian analogues of these 
conformal invariances. 

Although the coupling function Lo{(f)) is unconstrained by the structure of 
scalar-tensor gravity theories, its choice determines the form of the cosmo- 
logical models in the theory and the form of the weak-field limit. Typically, 
the weak-field solar-system predictions of scalar tensor gravity theories have 
the following relationship to those of general relativity (GR) 

a;(0) weak field result (GR result) x 



(48) 
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Specifically, the perihelion precession of Mercury, Ai), and the time- 
variation of G is predicted at second-order to be (Nordvedt 1970, Wagoner 
1970, Will 1993) 



per 100 yr (49) 



G _ /3 + 2cu\ 
G ~ ~ V4 + 3uj) 

where Go is the present measured value of G{t). There is even a special 
theory for which the term in [...] brackets vanishes in (50) and G{t) = Go is 
constant to this weak-field order (Barker 1978). In Brans-Dickc theory we 
have G oc 0^^. The general relativistic limit of an uj{(p) scalar-tensor theory 
is obtained (if it exists) by taking the two hmits 

uj' 

a; — > oo and ^ — > 0. (51) 

The observational hmits on uj require a; > 500 but the limit on u' is weak, 
Lo' < 0(1), (see Will 1993). Similar hmits are obtained from the binary pulsar 
(Damour, Gibbons & Tayler 1988), but are more model dependent. 



Ai? ~ 43" X 



1 - 



12 + Qlo 



4 + 



UJ 



(3 + 2ujy 



G{t) 2uj'{(f)) 



+ 



Go (3 + 2cu)2 







(50) 



4.1 Scalar- tensor Friedmann cosmologies 

In order to examine some of the counterparts to the Newtonian solutions 
discussed above we give the field equations for an isotropic and homogeneous 
Friedmann universe with scale factor a{t), curvature parameter k, and Hubble 
rate H = a/a, 

2 _ 8ttp Hip Lo{4))'(j)^ k 

+ 3i/0 + ^ = ^^^^^i^ (53) 
^ ^ 2a; + 3 2a; + 3 ^ ' 

p oc a-^'^. (54) 

We will now focus attention upon solutions of these equations with A; = 
in the case of Brans-Dicke theory (a; constant). Scalar-tensor theories with 
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varying G differ from general relativity in that they admit vacuum solutions 
when A; = (O'Hanlon & Tupper 1970) 

G oc (/.-^ocr'^/M (55) 
a{t) oc iV3(i+d) 

They also possess a class of special power-law solutions for perfect-fluid uni- 
verses, (Nariai 1969), 

a oc oc r'^ (56) 

a{t) oc 

P = il- 1)P 



where 

^ _ 2 + M2-T) (5 

4 + 3u;7(2-7) ' ^ ^ 

S . ^"7^^) ^ . (58) 
4 + 3a;7(2 - 7) ^ ^ 

The general solutions can be found for all 7 but are rather cumbersome and 
opaque; exact solutions for k = have been found by Gurevich, Finkelstein 
and Ruban (1973) and for all k by Barrow (1993a); a phase plane analysis 
has been performed by Kolitch and Eardley (1995) which includes the k ^ 
models. Their general properties are as follows. As i — > they approach the 
vacuum solutions (92), while as i — > 00 they approach the matter-dominated 
solutions (93)-(95). As c<j — > 00 these matter-dominated solutions approach 
the general relativistic results, a(t) oc t'^^^'^, G constant. There is a smooth 
transition between these simple early and late time behaviours. Thus the 
power-law matter-dominated solutions (93)- (95) are unstable as t — > 0. The 
general solutions are dominated by the Brans-Dicke scalar field. A similar 
early-time behaviour occurs in the cj(0) theories although the form of the 
early vacuum-dominated phase depends on the detailed functional form of 
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uj{(f)), (see Barrow 1993a, 1993b, Barrow and Mimoso 1994, Damour and 
Nordvedt 1993, Serna and Alimi 1996, and Barrow and Parsons 1996 for 
details) . 

4.2 New Methods of Solution 
4.2.1 Vacuum and radiation models 

The general solutions to Eqs. (52) - (54) contain four arbitrary integration 
constants, one more than their GR counterparts, the extra degree of freedom 
being attached to the value of 0. When the energy-momentum tensor is trace- 
free there exists a conformal equivalence between the theory and GR, the 
right-hand side of Eq. (53) vanishes and = is always a particular solution, 
corresponding to a special choice of the additional constant possessed by 
the model over GR. Consequently, the exact general solution of Einstein's 
equations when Tab is trace-free is also a particular solution to Eqs. (52)- 
(54) with (p, and hence a;(0), constant. 

It will seldom be the case that the particular solution obtained in this way 
will form the general solution for that particular choice oi u;{(j)) . Usually, 
however, it will be the late or early time attractor of the general solution. For 
example, in the case of Brans-Dicke theory the special GR solution is the late- 
time attractor for fiat and open universes but not the early-time attractor. 
However, a method has been developed for integrating the field equations 
for models containing trace- free matter (Barrow 1993a). The procedure is as 
follows. 

Eq. (54) integrates immediately to yield 



where F > is a constant of integration; F = describes vacuum mod- 
els. Making the choice 7 = 4/3, corresponding to blackbody radiation, and 
introducing the conformal time co-ordinate, rj, defined by 



8np = 3ra-^^ , 



(59) 



adrj — dt , 



(60) 



Eq. (53) becomes 




(61) 
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where subscript t] denotes a derivative with respect to conformal time. This 
integrates exactly to give, 

= 3^/2^(2cj(0) + 3)^^/2 ; A const. (62) 

We now employ the variable (suggested by the conformal invariance) used 
by Lorenz-Petzold (1984) to study Brans-Dicke models, 

?/ = (63) 

to re- write the scalar-tensor version of the Friedmann equation, Eq. (52), as 

{y,f = -Aky' + 4ry + ^ (0^)' a^(2o;(0) + 3) . (64) 

Dividing Eq. (|6^) by Eq. (|63D, and using Eq. (|62D, we obtain the coupled 
pair of differential equations 

^ = 3'/'Ay-\2u;{<P) + 3r'/\ (65) 

{y,f = -Aky^ + 4ry + A\ (66) 

We may now obtain the general solution for a particular choice of uj{(f)), given 
k. Integrating Eq. ( |65|) yields y{rj) which, in conjunction with co'(0), implies 
0(77) and, without further integration, 0(77) from Eq. (^). If Eq. ( |BUD can 
be integrated and inverted we may compute 0(t) and a(t), so completing the 
solution. The vacuum models are obtained by setting F = 0. 

4.2.2 General perfect-fluid cosmologies 

When T,the trace of the energy-momentum tensor is non-vanishing, the sit- 
uation is substantially more complicated. In this instance, = is no longer 
a particular solution of the field equations, forcing us to resort to more elabo- 
rate methods to obtain solutions. Barrow and Mimoso(1994) have done this, 
for the k = models, by generalising the method of Gurevich et al. (1973) 
for BD models to the case of non-constant a;(0). We outline this procedure. 

Introducing the new time co-ordinate ^, and the two new variables x and 
V such that 

dt = a'^^-^^ di , (67) 
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X = 



V = 



a 



3(2-7), 



(68) 
(69) 



and confining attention to the k = models, Eqs. (52)-(54) transform to 



2 \2 /2cj + 3 



and 



= 3r [(2 - 7)cu + 1] + 



3 / 2 



-a; + f c<j^ 



(70) 
(71) 

(72) 



2 V3(2cu + 3) 

wliere subscript-^ represents a derivative witli respect to ^-time. Eqs. (^ 
and (|72D integrate easily to yield 



r(4-37) (e-^i 

3 



X 



V2uj + 3 (^C + r(2 - 7) V2a; + 3 d^^ 



(73) 
(74) 



C is an integration constant and fixes the origin of ^-time. Noting the 
relation 

3 ' ^ ^ ' V i^^i _ f^V ^ (75) 



a(j) ''^'^ (f)^ ^'^^^ a _ ^ 
and differentiating y, with respect to ^, yields 



+ 



where a new function /(O? is defined by 



2r(4 - 37) 

Solving Eq. (|76|), we have the solution 



(76) 



C + r(2-7) / ^2cu(0) + 3rfe 



In 



(77) 
(78) 
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with (yf(^) simply related to /(^) by 



D. 



(79) 



where D is a constant of integration. Eq. (|7^) immediately reveals a simple 
formula for the scale-factor: 




2-7 



ao constant . 



(80) 



Finally, the scalar-tensor coupling function uj{(j)) is given as a function of / 
by 

^^m)) + 3 = 4^ r. ^{21 .1 , (81) 



where /o is another arbitrary constant. 

An important benchmark is provided by the behaviour of the BD theory, 
where uj{(j)) = uq = constant. In this case, the generating function, /(^), is 
given by a quadratic in ^: 



/BD(e) 



3(2-7) 
21(4 - 37) 



V2uJo + 3 



(82) 



we see that /bd oc ^ 
(23(7-i)g^^_ If we choose 



Hence, in general (C7 7^ 7^ F) when 7 7^ 4/3, 2 
as ^ ^ cxD and /bd oc ,^ as ^ ^0, where dt oc 
C = then /bd oc F^^ as ^ — » 0. The choice C = restricts the solution 
to the special 'matter-dominated' solutions (termed 'Machian' by Dicke, see 
Weinberg (1972)) which were first found for all perfect-fluids by Nariai (1969). 
If C 7^ then the early-time behaviour is dominated by the dynamics of the 
0- field; such solutions are termed '0-dominated' (or 'non-Machian' by Dicke). 
Therefore if we choose a generating function g{C,) that grows slower than 
as ^ — s> 00 it will produce a theory that approaches BD at late times 
{(f) ^constant, uj{4>) constant), whilst if g{^) decreases slower than ^ as 
^ — > then the theory will approach the behaviour of 0-dominated BD 
theory at early-times. This means that we will find new (non-BD) late- 
time behaviours by studying generating functions which increase faster than 
g{^) = as ^ ^ 00 and new (non-BD) early-time behaviour by picking 
generating functions which decrease slower than g{C,) = C, as ^ ^ or ,^ — >■ 
^min (if there is no zero of C, at the minimum of a{t)). 
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4.3 Comparing Newtonian and Relativistic Cosmolo- 
gies 

It is now possible to consider some of the similarities and differences that 
exist between the Newtonian cosmologies with varying G and their curved 

space-time counterparts. We recall that the special Newtonian power-law 
solutions (24) have the form G oc t""" with scale factor r oc t(2-n)/37_ ^j^jg jg 
identical to the Brans-Dicke solution when p = and 

n = — - — . (83) 

A + 3uj ^ ^ 

In fact, for any choice of a;(0) dust universes have 

(j)a^ (84) 
which corresponds to Newtonian solutions with G oc t^" and scale factor 

The vacuum solutions have a slightly different structure. If we identify 
then we have 

a{t) oc tV3(l+c^) 0, i(l-n)/3_ (gg) 

However, there is no general correspondence for other equations of state. 
Most notably, Brans-Dicke radiation-dominated universes (7 = 4/3) are the 
same as in general relativity, G = constant and a{t) oc t^/^, and differ from 
the radiation solutions. In general, the Newtonian solutions are the same as 
the Brans-Dicke matter-dominated solutions only if we make the choice 

4 + 3cJ7(2 - 7) ^ ^ 

Also, as can be seen from the right-hand side of the (j){t) evolution equa- 
tion, (90), the effective sign of the gravitational coupling changes sign from 
negative to positive when 7 becomes greater than 4/3. 

The Newtonian inflationary solutions with 7 = 0, given in cqn. (29), do 
not have direct counterparts in Brans-Dicke gravity theories. However, scalar 
tensor theories with 
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2a; + 3 oc t/)'' (88) 

have solutions of similar form as i — > oo with (Barrow and Mimoso 1994, 
Barrow and Parsons 1996) 

G oc oc t-2/{2^+i) (89) 
and the scale factor evolves as 

ait) oc t^f^-W2h+i) |^^2V(2?.+i) I _ (90) 

Thus, with 2/i + 1 > 0, we have a; — >• oo and uj'/uj^ — >• as t — >• oo 
and general relativity is approached in the weak-field limit. These solutions 
are not of identical form to the Newtonian solutions with G oc and 
n = 2/(2/1 + 1). 

5 Gravitational memory 

The process of primordial black hole formation in a cosmological model with 
varying G creates an interesting problem. We know (Hawking 1972) that 
black holes in scalar-tensor gravity theories are identical to those occurring 
in general relativity. Suppose, for simplicity, that a Schwarzschild black hole 
forms in the very early universe at a time tj when the gravitational coupling 
G{tf) differs from the value G{to) that we observe at the present cosmic 
time, to- This black hole will have an horizon size Rf — 2G{tf)M ^ tf 
when it forms. We now ask what happens to this black hole during the 
subsequent evolution of the universe as the value of G changes with time in 
the background imiversc (Barrow 1992, Barrow and Carr 1996). 

There appear to be two alternatives. The value of G{t) on the scale of 
the even horizon could change at the same rate as that in the background 
universe. But this would imply that the black hole was changing with time. 
Thus it could not be one of the black holes defined in general relativity, as 
required. Moreover, if G < 0,the horizon area (Ahor) would decrease with 
time and so would the associated entropy which is given by (Kang 1996) 

r, _ Ahor X 
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where (pi^^j. = G^J^ determines the value of G on the horizon. Ahcrnatively, 
there might be a process of 'gravitational memory' (Barrow 1992) wherein 
the scalar field determining G remains constant on the scale of the black 
hole horizon whilst changing in the cosmological background. That is at 
any moment of cosmic time there would be a space variation in G. This 
has dramatic implications for the Hawking evaporation of primordial black 
holes. For the lifetime and temperature of an evaporating black hole will be 
determined by the value of G{tf ) at the time when it formed rather than by 
the value G(to) we observe in the universe today: the black hole 'remembers' 
the value of G at the time of its formation. Hence, its Hawking lifetime will 
be ~ GjM^ and its temperature T^h ~ G^^M"^. 

Black holes which explode today are those whose Hawking lifetime is 
equal to the present age of the universe. This fixes their masses to be 



Clearly, if there has been significant change in the value of G since 
~ lO^^^s then the physical characteristics of exploding black holes can 
be very different from those predicted under the assumption that G does 
not change (Hawking 1974). Quite modest amounts of time evolution at 
unobservably early times can shift the spectral range in which we would see 
the evaporation products out of the gamma-ray band. This means that we 
should be looking for the evidence of black hole evaporation in other parts 
of the electromagnetic spectrum. A more detailed study of the observational 
evidence, seen is this light, is given by Barrow and Carr (1996). 



In recent years there has been a good deal of speculation about mechanisms 
for explaining the values of the fundamental constants. At one time it was 
widely believed that some ultimate Theory of Everything would eventually 



2 




and their temperature when the explode is therefore given by 




6 Origins of the Values of Constants 
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tell us that the constants could have one, and only one, set of logically self- 
consistent values. Such a simple scenario now seems less and less likely (Bar- 
row 1991). There are so many sources of randomness in the process which 
endow the fundamental constants with their low-energy values, and the pa- 
rameters likely to be fixed by a Theory of Everything are so far removed from 
our three-dimensional physical constants, that many new possibilities must 
be taken seriously. The non-uniqueness of the ground state of any Theory 
of Everything would mean that fundamental constants could take on many 
self-consistent sets of values. We would have to use anthropic constraints 
in order to understand those that we observe. This creates new interpreta- 
tional problems. Our underlying Theory of Everything would have quantum 
gravitational characteristics and its predictions about constants would have 
a probabilistic form. Although, formally, there would be a most probable 
value for the low-energy measurement of a quantity like the fine structure 
constant, such a value might be irrelevant for observational purposes (Bar- 
row 1994). We would only be interested in the range of values for which the 
evolution of complexity, in the form that we call 'life', is possible. This may 
well confine us to a subset of values which, a priori, is extremely improbable. 
This shows that in order to make a correct comparison of the probabilistic 
predictions of such a theory with observation we would need to know every 
dependence of processes which can lead to the evolution of complexity on the 
values of the constants of Nature. 

The fact that the observed values of many of the constants of Nature fall 
within a very narrow range for which life appears to be possible has elicited 
a variety of interpretations: 

(i) Good luck: the constants are what they are and could be no other 
way. The range which allows intelligent observers to evolve and persist is 
narrow and we are very lucky that our universe falls within that range. No 
matter how improbable this sate of affairs we could observe it to be no other 
way. 

(ii) Life is inevitable: we have been misled by our limited knowledge of 
complexity into thinking that life is restricted to universes spanned by a very 
narrow range of values for the constants. In fact, life may be a widespread 
inevitability in the phase space of all possible values of the constants. Even 
our own form of carbon-based life may exist in other novel forms which ex- 
ploit the possibilities provided by the recently discovered fuUerene chemistry, 
complexity of the sort that we call life may also exist in quite different forms 
to those we are accustomed to: for example, existing in velocity space rather 
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than in position space. 

(iii) All possibilities exist: whether through the actualisation of quantum- 
mechanical many-worlds, the realisation of all logically consistent Theories 
of Everything, or some elaboration of the self-reproducing universe scenar- 
ios, every possible permutation of the values of the constants exists in some 
universe. We live in one of the subset which allows life to exist. It is also 
possible that the ensemble of possibilities is played out in a single infinite 
Universe and we inhabit one of the life-supporting parts of it. 

(iv) Cosmic fine tuning: some physical process brings about approach 
to a particular set of values for the constants over long periods of time, 
perhaps through many cycles of cosmic evolution. The attracting set may be 
predictable in certain respects. 

The last of these four possibilities has attracted some interest recently. 
Harrison (1995) made the amusing suggestion that the fine tuning of the 
constants may be the end result of successive intelligent interventions by 
beings able to create universes in the laboratory (something discussed in 
the literature even by ourselves! Farhi and Guth (1987)). Aware that cer- 
tain combinations of the values of fundamental constants raise the probabil- 
ity of life evolving, and able to engineer these values at inception, succes- 
sive generations would tend to find themselves inhabiting universes in which 
life-supporting combinations obtained to high precision. Although Harrison 
refers to his as a 'natural selection' of universes, it is more akin to artificial 
selection, or forced breeding. 

Linde (1990) has proposed generalisations of the self-reproducing eternal 
infiationary universe in which the values of the fundamental constants change 
from generation to generation. Although unobservable, this scheme has the 
merit of being a by-product of the standard chaotic infiationary universe 
scenario. 

A third scenario of this sort, which has attracted a surprising amount of 
attention is that proposed by Smolin (1992) who suggested that a bounce, 
or quantum tunnelling, occurs at all final black-hole-collapse singularities 
which transforms them into initial singularities for new expanding universes. 
During this process the constants of Nature undergo small random changes. 
It is expected therefore that selection pressure will act so as to maximise the 
black holes produced in universes as time goes on (no weighting of the volume 
taking part in this reproduction process is introduced though, as is the case 
in the self- reproducing inflationary universes). Thus, if our universe is the 
result of the action of this selection process over many cycles of collapse 
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and re-expansion, in which the constants have lost memory of any initial 
conditions they may have had, then Smolin argues that we would expect 
to be near a local maximum in the black-hole production. Hence, small 
changes in the constants of Nature should in general take us downhill from 
this local maximum and always reduce the amount of black hole production. 
By conducting such thought-experiments the general consistency of the idea 
can be tested. 

We make three remarks about this speculative scenario. First, it is not 
clear that is as sharply predictive and testable as claimed. We should only ex- 
pect to find ourselves residing near a local maximum in the space of constants 
if that maximum also provides conditions which permit living observers to 
exist. If those condition are unusual then we might have to exist in one of 
the improbable universes far from the local maxima. We can only deter- 
mine if this is the case by having a complete understanding of the necessary 
and sufficient conditions for living complexity to exist. Second, putting this 
objection to one side, there may well be small changes in the values of the con- 
stants which significantly increase the production of black holes. For exam- 
ple a small {70KeV) strengthening of the strong interaction would bind the 
dineutron and the diproton (helium-2), so providing a direct H + H ^ He^ 
channel for nuclear burning, massive stars would run through their evolu- 
tion very rapidly and end as black holes far sooner and with higher likelihood 
than at present (see Dyson 1971, Barrow 1987). Third, we might ask why 
there should be any local maxima at all for variations in certain constants. 
Variation would proceed to states of higher gravitational entropy by always 
increasing the value of Sbh oc GM^, and this would be effected by a random 
walk upwards through (over long time averaged) increasing values of G. 

Finally, we should note that any scheme which relies upon random changes 
in the constants of Nature occurring at the endpoint of gravitational collapse 
must beware of the consequences of changes which prevent future collapses 
from occurring. A specific example is seen in the case of closed universes 
oscillating under the requirement that their total entropy increase from cycle 
to cycle. There, one finds that any positive cosmological constant (no matter 
how small in value) which remains constant (or falls slowly enough on aver- 
age) from cycle to cycle ultimately stops the sequence of growing oscillations 
and leaves the Universe in a state of indefinite expansion which asymptotes 
towards the de Sitter state (Barrow and Dabrowski 1995). In Smohn's sce- 
nario one might consider that if the curvature of space or the value of the 
cosmological constant, or the magnitude of vacuum stresses associated with 
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scalar fields which violate the strong energy condition, were to change at the 
collapse event in ways that prevented future collapse of some or all of the 
Universe, then gradually the fraction of the Universe which could gravitation- 
ally collapse and evolve the values of its constants by random reprocessing 
would shrink asymptotically to zero. Evolution would cease. This Universe 
would have 'died'. 

7 Simultaneous Variations of Many Constants 

The subject of varying constants is of particular current interest because 
of the new possibilities opened up by the structure of unified theories, like 
string theory and M-theory, which lead us to expect that additional compact 
dimensions of space may exist. Although these theories do not require tradi- 
tional constants to vary, they allow a rigorous description of any variations to 
be provided: one which does not merely 'write in' the variation of constants 
into formulae derived under the assumption that they do not vary. This self- 
consistency is possible because of the presence of extra dimensions of space 
in these theories. The 'constants' seen in a three-dimensional subspace of 
the theory will vary at the same rate as any change occurring in the extra 
compact dimensions. In this way, consistent simultaneous variations of dif- 
ferent constants can be described and searches for varying constants provide 
a possible observational handle on the question of whether extra dimensions 
exist (Marciano, 1984, Barrow 1987, Damour & Polyakov 1994). 

Prior to the advent of theories of this sort, only the time variation of 
the gravitational constant could be consistently described using scalar-tensor 
gravity theories, of which the Brans-Dicke theory is the simplest example. 
The modelling of variations in other 'constants' was invariably carried out by 
assuming that the time variation of a constant quantity, like the fine structure 
constant, could just be written into the usual formulae that hold when it is 
constant. One way of avoiding this situation is to exploit the invariance 
properties of the non-relativistic Schrodinger equation for atomic structure, 
which allow it to be written in dimensionless form when atomic ('Hartree') 
units are chosen. It can be shown (Barrow and Tipler 1986) that any solution 
with an energy eigenvalue E, arising when the fine structure constant is a 
and the electron mass is me, must be related to a solution defined by a E, a', 
and rrig by the relation 
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where c is the velocity of hght. 

The possibihty of hnked variations in low-energy constants as a result 
of high-energy unification schemes has the added attraction of providing a 

more powerful means of testing those theories (Marciano 1984, Kolb, Perry, 
& Walker, 1986, Barrow 1987, Dixit & Sher 1988, Campbell & Olive 1995). 

Higher-dimensional theories typically give rise to relationships of the fol- 
lowing sort 

ai{m^) = AiGml = BiX^{ipi/R)'' -,11, k constants (92) 

tt"^ ^Cij[ln(m*/mj) + 9{iJ, - rrij) \n{mj/ jj)] + Aj 

where ai{..) are the three gauge couplings evaluated at the corresponding 
mass scale; ^ is an arbitrary reference mass scale, is a characteristic 
mass scale defining the theory (for example, the string scale in a heterotic 
string theory); A is some dimensionless string couphng; l^i — G'^l"^ is the 
Planck length, and i? is a characteristic mean radius of the compact extra- 
dimensional manifold; are numbers defined by the particular theories 
and the constants Ai and Bi depend upon the topology of the additional 
( > 3) dimensions. The sum is over j = leptons, quarks, gluons, ^ Z 
and apphes at energies above ~ \GeV (Marciano 1984). The term Aj 
corresponds to some collection of string threshold corrections that arise in 
particular string theories or an over-arching M theory (Antoniadis & Quiros 
1996). They contain geometrical and topological factors which are specified 
by the choice of theory. By differentiating these two expressions with respect 
to time (or space), it is possible to determine the range of self-consistent 
variations that are allowed. In general, for a wide range of super-symmetric 
unified theories, the time variation of different low-energy constants will be 
linked by a relationship of the form (where we consider /3 to denote the time 
derivative of /5 etc.) 

P G ^ af ^ rrij A 

where /3 = me/rripj. (Drinkwater et al 1997). It is natural to expect that all 
the terms involving time derivatives of 'constants' will appear in this relation 
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unless the constant 5 prefactors vanish because of supersymmetry or some 
other special symmetry of the underlying theory. This relation shows that, 
since we might expect all terms to be of similar order (although there may 
be vanishing constant 6 prefactors in particular theories), we might expect 
variations in the Newtonian gravitational 'constant', G/G, to be of order 

8 Varying alpha — New Observational Limits 

Quasar absorption systems present ideal laboratories in which to search for 
any temporal or spatial variation in the assumed fundamental constants of 
Nature. Such ideas date back to the 1930s, with the first constraints from 
spectroscopy of QSO absorption systems arising in the 1960s. An historical 
summary of the various propositions is given in Varshalovich & Potekhin 
(1995) and further discussion of their theoretical consequences is given in 
Barrow & Tipler (1986). 

Recently, we have considered the bounds that can be placed on the varia- 
tion of the fine structure constant and proton g factor from radio observations 
of atomic and molecular transitions in high redshift quasars (Drinkwater et 
al 1997). To do this we exploited the recent dramatic increase in quality of 
spectroscopic molecular absorption at radio frequencies, of gas clouds at in- 
termediate redshift, seen against background radio-loud quasars. Elsewhere, 
we will consider the implications of simultaneous variations of several 'con- 
stants' and show how these observational limits can be used to constrain 
a class of inflationary universe theories in which small fluctuations in the 
flne-structure constant are predicted to occur. 

The rotational transition frequencies of diatomic molecules such as CO are 
proportional to k/ {Ma?) where M is the reduced mass and a — Ti? / {vriee^) is 
the Bohr radius. The 21 cm hyperflne transition in hydrogen has a frequency 
proportional to ^^^^/{ha^), where fi^ = gpch / {Ampc) , Qp is the proton g- 
factor and ii^ = eh/{2mec). Consequently (assuming rUp/M is constant) 
the ratio of a hyperflne frequency to a molecular rotational frequency is 
proportional to Qpa^ where a — e^/ (he) is the flne structure constant. Any 
variation in y = gpO^ would therefore be observed as a difference in the 
apparent redshifts: ^z/{l + z) ^ Ay/y. Redshifted molecular emission is 
hard to detect but absorption can be detected to quite high redshifts (see 
review by Combes & Wiklind, 1996). Recent measurements of molecular 
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absorption in some radio sources corresponding to known HI 21 cm absorption 
systems give us the necessary combination to measure this ratio at different 
epochs. 

Common molecular and HI 21cm absorptions in the radio source PKS 
1413+135 have previously been studied by Varshalovich & Potekhin (1996). 
They reported a difference in the redshifts of the CO molecular and HI 21 cm 
atomic absorptions which they interpreted as a mass change of AM/M = 
(—4 ± 6) X 10~^ but as we show above this comparison actually constrains 
gipO;^, not mass. Furthermore they used overestimates of both the value 
and error. They used the Wiklind & Combes (1994) measurement which 
had the CO line offset from the HI velocity by —11 kms~^; a corrected CO 
measurement (Combes & Wiklind, 1996) shows there is no measurable offset. 
Furthermore, Varshalovich & Potekhin (1996) used the width of the HI line 
for the measurement uncertainty. Even allowing for systematic errors the true 
uncertainty is at least a factor of 10 smaller so these data in fact estabhsh a 
limit of order 10~^ or better. This potential for improved limits has prompted 
the present investigation: previous upper limits on change in a are of order 
Aa/a pa 10""^ (Cowie & Songaha, 1995; Varshalovich, Panchuk & Ivanchik, 
1996). 

8.1 Comparison of HI and molecular systems 

Our new more accurate rcdshift estimates (Drinkwater et al 1997) give a 
molecular redshift of 0.684680±0.000006 and an a 21cm redshift of 0.684684± 
0.000006 for the source 0218 + 357 and a molecular redshift of 0.246710 ± 
0.000005 and a 21cm redshift of 0.246710±0.000004 for the source 1413+135. 
We can therefore combine the uncertainties in quadrature to give 1-sigma 
upper limits on the redshift differences. These give |A^/1 + ^| < 5 x 10~^ 
(1.5 kms^^) for both sources. 

We must still consider the possibility that the molecular and atomic ab- 
sorption arises in different gas clouds along the hue of sight. This could 
explain any observed difference. However there is no measurable difference 
between the two velocities in our data, so we arc probably detecting the same 
gas. The alternative would be that there was a change in the frequencies but 
that in both cases it was exactly balanced by the random relative velocity of 
the two gas clouds observed. We consider this very unlikely because of the 
small lKms~^ dispersion within single clouds. 

We can now use the limits to Az/1 + z with the relationship to derive 
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1-sigma limits on any change in y = QpO^: \/\y/y\ < 5 x 10~^ at both z = 0.25 
and z = 0.68. These are significantly lower than the previous best limit of 
1 X 10""^ by Varshalovich & Potekhin (1996) (it was quoted as a limit on 
nucleon mass, but it actually refers to QpO^). 

As there are no theoretical grounds to expect that the changes in Qp and 

are inversely proportional, we obtain independent rate-of-change limits of 
\9pl9p\ < 2 X 10~^^y"-^ and \a/a\ < 1 x 10"-^^ y~^ at 2; = 0.25 and \gp/gp\ < 
IxlO^^^y-i and \a/a\ < 5x10-^^-^ at z = 0.68 (for iJo = 75kms-^Mpc-^ 
and qo = assumed here). These new limits are stronger than the previous 
1 sigma limit of |q;/q;| < 8 x lO^^^y"-^ at z ^ 3 (Varshalovich et al. 1996). 

The most stringent laboratory bound on the time variation of a comes 
from a comparison of hyperfine transitions in Hydrogen and Mercury atoms 
(Prestage et al. 1995), \a/a\ < 3.7 x and is significantly weaker 

than our astronomical limit. The other strong terrestrial limit that we have 
on time variation in a comes from the analysis of the Oklo natural reactor 
at the present site of an open-pit Uranium mine in Gabon, West Africa. 
A distinctive thermal neutron capture resonance must have been in place 
1.8 billion years ago when a combination of fortuitous geological conditions 
enriched the subterranean Uranium-235 and water concentrations to levels 
that enabled spontaneous nuclear chain reactions to occur (Maurette 1972). 
Shlyakhter (1976, 1983) used this evidence to conclude that the neutron 
resonance could not have shifted from its present specification by more than 
5 X over the last 1.8 billion years and, assuming a simple model for 

the dependence of this energy level on coupling constants like a, derived a 
limit in the range of \a/a\ < (0.5-1.0) x 10~^'^y~^. The chain of reasoning 
leading to this very strong bound is long, and involves many assumptions 
about the local conditions at the time when the natural reactor ran, together 
with modelling of the effects of any variations in electromagnetic, weak, and 
strong couplings. Recently, Damour & Dyson (1997) have provided a detailed 
reanalysis in order to place this limit on a more secure foundation. They 
weaken Shlyakhter's limits shghtly but give a 95% confidence hmit of —6.7 x 
10~^^y~^ < a/a < 5.0 x 10~^^y~^. However, if there exist simultaneous 
variations in the electron-proton mass ratio this limit can be weakened. 

These limits provide stronger limits on the time variation of a than the 
astronomical limits; however, the astronomical hmits have the distinct advan- 
tage of resting upon a very short chain of theoretical deduction and are more 
closely linked to repeatable precision measurements of a simple environment. 
The Oklo environment is sufficiently complex for significant uncertainties to 
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remain. 

Unlike the Oklo limits, the astronomical limits also allow us to derive up- 
per limits on any spatial variation in a. Spatial variation is expected from the 
theoretical result that the values of the constants would depend on local con- 
ditions and that they would therefore vary in both time and space (Damour 
& Polyakov 1994). The two sources for which we derived hmits, 0218-1-357 
and 1413+135, are separated by 131 degrees on the sky, so together with the 
terrestrial result, we find the same values of a to within |Aa/Q;| < 3 x 10~^ 
in three distinct regions of the universe separated by comoving separations 
up to 3000 Mpc. Limits on spatial variation of gpO^me/mp were previously 
discussed by Pagel (1977, 1983) and Tubbs & Wolfe (1980). We have im- 
proved on their limits by some 2 orders of magnitude but as our sources are 
at lower redshift, they are not causally disjoint from each other. 

The high-redshift measurements are now approaching the best terrestrial 
measurements based on the Oklo data. These could be further improved 
by a factor of 2-5 with additional observations that would not be difficult 
to perform such as fitting the atomic and molecular data simultaneously, 
remeasuring the HI absorptions at higher spectral resolution. 

8.2 Inflation 

Inflation is something of a two-edged sword when it comes to discussing 
variations in constants. On the one hand there are potentials with multiple 
vacuum states which allow different parts of the universe to find themselves 
inheriting different suites of fundamental constants, with quite different val- 
ues. On the other hand, if, as inflation leads us to expect, the whole of 
our observable universe is contained within the inflated image of a single 
causally connected region, then we should expect fundamental constants to 
reflect that single origin and to display spatial uniformity to very high pre- 
cision. The key question is what precisely is that precision? A bench mark 
for the amplitude of possible variations is provided by the amplitude of tem- 
perature fluctuations in the microwave background, AT/T ~ 10~^. We 
would expect fluctuations in the fine structure constant created at the end 
of inflation to have an almost constant curvature spectrum (because of the 
time-translation invariance of almost de Sitter inflation) with an amplitude 
below that of 10~^. An interesting feature of the new astronomical observa- 
tions described above is that, for the first time, they take the observational 
limits on spatial variations in a ( |Aq!/q;| < 3 x 10~^) into that regime where 
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they may be constraining the underlying theories more strongly than are the 
COBE observations. 
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